SOME NONLINEAR FUNCTIONS OF BERNOULLI AND EULER UMBRJE 



CHRISTOPHE VIGNAT, UNIVERSITE D'ORSAY 



Abstract. In a recent paper [5], Yi-Ping Yu has given some interesting nonlinear moments of the 
Bernoulli umbra; the aim of this paper is to show the probabilistic counterpart of these results and 
to extend them to Bernoulli polynomials. 

b 
a . 

1. Introduction 

In a recent rich contribution, Yi-Ping Yu gives several nonlinear moments of the Bernoulli umbra 
05 defined by its generating function 

exp(zOS) = -A: -, |*| < 2vr. 

> . exo (z — 1 



exp {£) — 1 ' 

This umbra is related to the Bernoulli numbers as 

93" = B n ; 

for example 

Bo = 1; _Bi = — ; Bo = — ; B?. — 0; Ba = . 

0,1 2 6 30 

and all odd orders Bernoulli numbers except B± equal 0. 

Similarly, the Euler umbra € is defined by the generating function 

exp (z<£) — sech (z) 

We generalize here these umbrae and define the Bernoulli umbra 55 (a;) as 

(1.1) exp (2*8(3)) = -^- 

e z — 1 



u 



> 
o 



and the Euler umbra 2; (x) as 



exp (a;)) = - 



As a result 
and 



f 



<T (x) = ^ n (x) , 

respectively the Bernoulli and Euler polynomials of degree n. 

The aim of this paper is to compute some nonlinear functions of these umbrae as probabilistic 
nonlinear moments. In the following, we denote the expectation operator 

Eh {X) = J h (x) fx (x) dx 

where fx is the probability density function of the random variable X. We will use the following 
characterization of the Bernoulli and Euler umbras. 

Theorem 1. The Bernoulli umbra 05 (x) satisfies, for all admissible function h, 

h(&(x)) = Eh(x- \+iLb 
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where the random variable Lb follows a logistic distribution, with density 

Il b (x) = -sech (nx) , x G M. 
Accordingly, the Euler umbra (£ (x) satisfies, for all admissible function h, 

h (<£ (x)) = Eh(x- -+lL E 



2 

where the random variable Le follows the hyperbolic secant distribution 

Il e i x ) = sech ( nx ) ■ 

Proof. Since 

exp (rf*8 (x)) = Ecxp ^it ^x — — + iLb 

by identification with , the random variable Lb has characteristic function 

t 

n itL B \ _ 2 



V ' sinh(|) 
But from 1.9.2] 

/ sech 2 (ax) cos (xt) dx = -^4-csch f — 
Jo 2a 2 \2a 

so that, with a = n, the density of Lb is 

Il b (x) = |sech 2 (nx) , 

which is a logistic density. 

Accordingly, the characteristic function of the random variable Le is 

't 



Ee lLEt = sech 

From [6, 1.9.1], 



+ oo 







7T / 7T \ 

sech (ax) cos (xt) dx — — sech ( — t 
so that, with a = 7r, the density of L E is 

Jl e ( x ) = sech(7rx) . 

Thus ttLq follows an hyperbolic secant distribution. □ 
As a consequence, the Bernoulli polynomials read 

(1.2) B n (x) = 58 (x) n = E (x - i + iL B ^J 



and the Bernoulli numbers 



B n = <8 n = <B (0)" = E ( -\ + iL B ) , n > 0. 



Similarly, the Euler polynomials read 

E n (x) = € (x) n = e(x-^+iL e 

and the Euler numbers n 

E n = 2 n z(^j =2 n E(iL E ) n . 

We note from [7J p. 471] that the random variable Lb can also be obtained as 

L H = — log = — log 

2tt 6 1 - U 2vr 6 E 2 
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where U is uniformly distributed on [—1, +1], E\ and Ei are independent with exponential distribution 
Je (x) = exp (—x) , x € [0, +oo[ and equality is in the sense of distributions. As for the random variable 
Le, from [TJ, it can be obtained as 

(1.3) Lo = - log |C| - - (log \Nx\ - log |JV 2 |) 

7T 7T 

where C is Cauchy distributed and N\ and N2 are two independent standard Gaussian random vari- 
ables. 



We compute 



which, by symmetry, is equal to 



2. THE MOMENT log 03 (x) 



logQS (x) = E\og[x - - + iL B 



ifilognar-i) / I =log 



1 , ( L 2 R 
-E log 1H 



but from O 2.6.30.2] 



» sinh cz 

Thus, by bisection of the angle 2irz, 



+°°log(l + 6z 2 ) d 2 1 c 

dz — n(b,c) = — log — j= — ip 



ny/b \ix\fb 



+°° log(l + &z 2 ) 1 f + °° \og(l + bz 2 ) 1 f +0 ° log (1 + bz 2 ) /cosh 2 TTZ 

^— 2 -dz = — / ^ ^ — = — / — ^ L ^ 1 ] dz 

, 1 sinh 2-kz 4 J sinh 7rz cosh 7rz 4 J cosh nz \ sinh 7rz 



so that 



+°° logfl + 6z 2 ) 

6V 2 ^-dz = n(h(b,Tr) -4h(b,2n)). 



00 



COSh 7TZ 



We deduce, with b = (x — |) 2 , 



7-2 



-Slog 1 



and, using the identity 



= f (M M)-4MM*)) = f0(ioi^-^ 
" bg 7E- 21og ^- li '(7E) +2 *(7? 



^ (22) = i v (2) + ^ u + ~ ) i< "A 2. 



we obtain after simplification 



Slog (iL B + x - ~J = log -L + Slog (1 + 6L|) = i, Q 
For x — 1, we recover the result by Y.-P. Yu, namely 



.x 

1 2 



n..g [ I + )=,-(!) = --. 
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3. THE MOMENT log <& (x) 

This moment can be obtained according to the same approach, namely, again with b = (a; — ^) 

log (S (x) = log -j= + X -E log (1 + bL%) 
where the latter expectation is now computed using [3, 2.6.30.1] as 

+°° log (1 + &z 2 ) r (f + i75) . 1 

^— — r-t-dz = 2 log — t (- - log — = 

o cosh(Trz) r(± + -kS\ 2Vb 

V 4 2VT 



so that 



log£(x) = log 2 



2VbJ 

r 2 (l + il^-^l) 



4. THE MOMENTS *B k (x) AND € k (x) 

By derivation of the preceding results, we deduce 

2T 1 (x) =e(x-±+ iL b ^ = A bg © (z) 

so that we have 

(ip'(x), x>\ 
OS" 1 (a) = I -ip' (-x + 1) , x<i 
[o x = i 

and we remark that Q^ 1 (a;) is not continuous in x = \. Since moreover for any integer k > 1 

•'• 2 + H = (fc- l)l^ E { X -2 +lLB 



we deduce 



V 2 y ^_ I _2_^W(_ a; + i) a;<i 

and in a particular case x = 1, since ^( fc )(l) = (-l) fe+1 fe!C(fc+l) ! 

»-*(!) =-eQ + iLb) =fc.C(fc + l). 



In the Euler case, we have 



z = A 



and € 1 (x) is not continuous in a; = ^. 
More generally, for any integer k > 1, 



l 
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5. THE MOMENT logsin 2 ^ 

This moment can be easily computed from the moment representation as follows 

7rQ3 7T / 1 \ / 7T TtLb\ ( 7T 7tLb 

logsm — = E log sin — I — + iLg I = E log sin I - i - I = E log sin I 7 + 4 



2 2 v 2 / °V 4 2/ °v 4 2 

1 „, . / 7T 7T_L S \ . / 7T 7rL B 

= -£/ logsm sin — — + «— — - 

2 S V 4 2 / \ 4 2 

and expanding the product of sines we obtain 

Tplog Q cos ) + ^cos(i7ri B )^ = -ilog2 + ^E log cosh (ttL b ) . 

But since 

^ = 2 l0g T=u> 

we deduce 

cosh (itLb) 



2y/U(l-U) 

so that, with E'log U = — 1, we deduce 

i? log cosh (ttLb) = — log2 + 1 



and the result 
follows. 

The Pochhammer symbol 

has been evaluated in [TJ p. 149] as 



• ttQS 1 , „ 
log sin — = - - log 2 



6. the Pochhammer (*8 (x)) 7 
r(*B + n + l) 



(23 + 1) 



r(«8 + i) 

11 1 

(» + !)« 



(n + 1) 

We use the "intuitive argument" suggested by Carlitz [2j to compute its polynomial version (23 (x)) 7 
as follows: a generating function of (03 (x)) is 

p(a;,t) = ^(B(a;)) n ^=f;exp(-(a:-i+tL^tog(l-i) 

= (1 -fT*^) £:exp(-« J L B log(l -*)) 
with the characteristic function for the logistic density 

u 

Eexp(iL B u) = , 2 , M , 
smh(|) 

so that 

This term is identified as the derivative 

d (i - ty^-v _ d ^ t"- 1 , _ _v^^Ar 

i ~ dx ^ nl ( ' T 1)n ~ ^ n! (X 1)n 

n=0 n=0 

with 

^-{x- 1)„ = (a: - 1)„ (V (x + n - 1) - ^ (a: - 1)) 
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so that the coefficient of ^ in ip (x, t) is 

(® (»))„ = (a: ~+r +i ( ^ (a; + n) - v ' {x i)) ■ 

We recover the result by Norlund by taking the limit — » 1 which is -frr. 

7. THE POCHHAMMER (<t(x)) n 

We use the same approach to compute the Pochhammer symbol of the Euler polynomial umbra; 
the generating function reads 

+°° £u / / 1 \ 

<p{x,t) = Y l (^(x)) n -=Eexp(-(x--+iL E Uog(l-t) 

= (l-ty( x ~^ Eexp {-iL E log (1- t)) 
with the characteristic function of the hyperbolic secant distribution 

Ee lLEt = sech ( - 
V2 

so that 



Eexp(iL E log(l -<)) = sech ( ilog(l-t) * v ' ' 



and 



so that 



1 - i 

A 2 



(i-*)- 1 (i-l) 



Z_j, n j 2 " ^— ^ fc! 

n=0 fe=0 



fe=0 
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